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Abstract. The Sc galaxy M33 has been mapped with ISOPHOT in the far-infrared, at 60, 100, and 170 yum. The spatial 
resolution of these FIR maps allows the separation of spiral arms and interarm regions and the isolation of a large number 
of star-forming regions. The spectral energy distribution in the FIR indicates a superposition of two components, a warm one 
originating from dust at ~ 45K, and a cold one, at ~ 16K. The warm component is concentrated towards the spiral arms 
and the star-forming regions, and is likely heated by the UV radiation from OB stars. The cold component is more smoothly 
distributed over the disk, and heated by the diffuse interstellar radiation. For the about 60 star-forming regions detected the 
HoyFIR flux ratio increases significantly with the distance from the galaxy center, probably due to decreasing extinction. An 
| anti-correlation of F Ha /F 60 with F 1J0 suggests the intrinsic extinction to be related to the cold dust surface brightness according 

, to Ay /S no ~ 0.03 magMJy -1 sr. For the total galaxy the star formation rate (SFR) derived from the FIR is in agreement with 

\^ ' that derived from the de-extincted Ho- emission. For individual star-forming regions, a consistency between SFRs derived from 

the optical and from the FIR requires only a fraction of the UV radiation to be absorbed locally. The individual star-forming 
CO . regions also show a local radio-FIR correlation. This local correlation is, however, due to quite different components than to 

' those that lead to the well-known global radio-FIR correlation for entire galaxies. 
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' 1. Introduction bulk of the FIR emission, as in the case of normal galaxies. 

& This cold component arises from grains which are in weaker 

• • The sensitivity of ISO (Kessler et al. 1996) and its spectral cov- radiation fiddg> ^ fa ^ Quter optically ^ regions of ^ 

£ erage extending to 200 ^m made it the first observatory capable galaxy diskfj Qr becaufje they ^ from radiation fey Qp _ 

^ . of measuring the bulk of starlight absorbed by dust in galaxies. ^ depth effects (Tuffs & popescu 2QQ2) 

5h ' Recent FIR studies of normal galaxies (e.g. Haas et al. 1998 for 

M 3 1 ; see also Tuffs & Popescu 2002 for a review), based on 

data taken with the ISOPHOT instrument (Lemke et al. 1996) r„, , . . e . , , , , 

v ' Though the existence of the warm and cold components 

onboard ISO, have proven the existence of a cold dust compo- , • , . t A e . .. .. , . ,■ 

' F F were unambiguously inferred from the statistical studies men- 

nent with a temperature below 20 K, in addition to the warm .• , , .. , . t - ftU 

F ' tioned above, the spatial distributions of these components can 

component near T ~ 50K. These components were originally , , . , , , • f , , , , . 

F F 6 1 only be studied by the analysis of resolved nearby galaxies, 

proposed by Chini et al. (1986), based on sub-mm observa- Such an analyfjis waf , done for ^ IS OPHOT map of M 31 by 

tions. Statistically significant evidence for the existence of a Haaf , gt aJ (199g) and Schmidtobreick et al . (2000) , who re . 

cold dust component was also inferred from the ISOPHOT solved seyeral dozens of bright ^ and separated their SEDs 

Serendipity Survey (Stickel et al. 2000), and from the analy ^ cQ j d and warm dust com p 0nen ts. While the FIR emission in 

sis of a complete volume- and luminosity limited sample of cu , .,. ,. ... •, 

F J F M 31, a Sab galaxy seen at high inclination, originates mainly 

late-type Virgo Cluster galaxies - the ISOPHOT Virgo Cluster from a rf ^ stmcmre rather than from ^ ^ M a 

Deep Sample (Tuffs et al. 2002, Popescu et al. 2002). Even c , „ .. . . , ■ , . e 

f f x ' Sc galaxy with well distinct spiral arms, and seen almost face 

ultra-luminous IR galaxies show evidence for a cold dust com- . , , . . . . , .. , .■ 

& on, seems to be an ideal target to study the spatial distribution 
ponent (Klaas et al. 2001), though this does not constitute the c , . , .. , • . , , . ^ . .. 

F v " 6 of dust and its heating sources in spiral galaxies. Due to its 

* Based on observations with the Infrared Space Observatory ISO. distance of only 830 kpc, its scale of 0.24 kpc arcmin- 1 is suffi- 
ISO is an ESA project with instruments funded by ESA Member cient to allow a separation of its constituents (nucleus, Hn re- 
States (especially the PI countries France, Germany, The Netherlands gions, spiral arms, and interarm regions) with ISOPHOT even 
and the UK) and with the participation of IS AS and NASA. at 170 fim. 
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We will make use of the higher spatial resolution and the 
extended spectral coverage (to longer FIR wavelengths) of 
ISOPHOT as compared with IRAS to disentangle the warm 
and cold dust components in M33. The spatial distributions of 
these dust components will be compared with the distribution 
of the optical emission. The FIR emission from individual HII 
regions will be correlated with the corresponding radio emis- 
sion. In this paper we will also study the relation between FIR 
emission and star formation rate, as well as the heating of the 
two dust emission components and their effects on extinction. 



2. Observations and data reduction 

2.1. FIR data 

M33 was mapped with ISOPHOT (Lemke et al. 1996) in the 
raster scan mode, AOT-P32, at three wavelengths, 60, 100, and 
170/im, using the CI 00 and C200 detector arrays. The raster 
step sizes in spacecraft coordinates were 69" x 92" for 60 and 
lOOyum and 92" x92" at 170yum. Combined with the detector 
geometry and the motion of the focal plane chopper, this yiel- 
ded a sky sampling of 15" x23" and 31" x92" respectively. 
Since the linear extension of M33 is too large to be coverable 
by one raster map with the desired resolution, two partial maps 
had to be done in each filter (North and South) and concate- 
nated. With the partial Northern and Southern P32 maps mo- 
saiced together according to their center positions and satellite 
roll angles, a total area of 48' x 32' was covered, including the 
entire disk of M33, and a surrounding background area to al- 
low for proper sky subtraction. The orientation of this map is 
roughly along the major axis of the galaxy. 

At 170/mi the central area of 29' x 26' size was also mapped 
in the AOT-P22 raster staring mode, which almost is free from 
detector transient effects due to the relatively long integrating 
time of 10 seconds per raster position. The P32 and P22 maps 
were merged to a common one by using equal weights for both 
maps in the bright and medium bright areas and zero-weight for 
the P22 map for areas where the surface brightness falls below 
13MJy/sr. 

The data reduction was done using the ISOPHOT 
Interactive Analysis software package PIA V7 (Gabriel et al. 
1998), together with the calibration data set V4.0 (Laureijs et 
al. 2000). A special data handling was required due to signal 
transient effects of the infrared detectors, which is described 
in sect. 2.2. The absolute flux calibration was performed with 
measurements of the internal thermal fine calibration source on 
board, before and after each map. We estimate the error of this 
flux calibration to be about +20%. An additional uncertainty 
is introduced by the detector transient effects in the P32 ob- 
serving mode, which generally leads to too low signals for the 
FIR knots, and which we have only corrected to a certain de- 
gree (see below). At 60 and 100 /im, where this effect is most 
critical, we estimate the total error to be +30%. At 170yum the 
transient effect is less severe, and the total error is of the order 
of ±20%. 



2.2. Handling of transient detector effects 

As mentioned above, a difficulty in the reduction of P32 data 
are signal transient effects of the infrared detectors, leading to 
small lags of the maxima positions for bright point-like sources 
along the scan direction, and to a lowering of the maxima and 
filling up of the minima. It also led to ghost images of bright 
sources at distances of the chopper throw in the scan direc- 
tion. Based on the principles described by Acosta-Pulido et al. 
(2000) we adopted a correction routine to reduce the transient 
effects: 

First, separate maps were constructed for each detector 
pixel, nine for the 60 and 100//m maps, which were observed 
with the CI 00 detector array, and four in the case of the 170/im 
map, which was observed with the C200 array. Transient ef- 
fect were reduced by adding to the signal in each pixel 30% of 
the difference to the signal of the preceeding pixel during the 
scan. This procedure led to a steepening of the edges of bright 
sources and improved the coincidence in the forth and back 
scans. It also led to a significant reduction of the ghosts from 
bright point like sources at the map pixels located in chopper 
throw separation (here, however, a factor of only 10% was ap- 
plied to the difference signal). 

Second, the separate pixel maps were scaled to a common 
level and superimposed. Due to a displacement of the raster 
legs by half of the array size and a chopper step of 15" the pixel 
size of the resulting image was 15"x23" (in the cases of 60 and 
100//m). Each of these pixels was seen by 3 detector pixels in 
the one raster leg in forth direction and in the adjacent raster 
leg in opposite direction by 3 other detector pixels. From the 
6 detector signals in each pixel the median signal was derived, 
and an iterative k-ct clipping was applied toin the separate pixel 
maps and to correct them. This also leads to a spatial resolution 
of 60" at 60//m, of 80" at lOOyum, and of 120" at 170/rni. It was 
found that each scan had a slight positive slope, recognizable 
by small offsets at the edges of the forth and back scans. This 
was corrected by a fine adjustment of the scan slopes until a 
smooth background was achieved. The background level (from 
foreground cirrus and zodiacal light) was then determined from 
the outermost map areas and subtracted. 

2.3. Optical observations 

For studying the correlation of optical emission with the de- 
tailed ISOPHOT maps, M 33 was observed in Ha with the focal 
reducer at the 2.2 m telescope on Calar Alto, Spain. The detec- 
tor was a 2kx2k SITe chip used in a rebinned mode, to provide 
a scale of 1'.'06 pixel -1 , thus covering a field of 15'diameter. 
The photometric conditions were good. The seeing of 275 was 
adequate for this study. For isolating the Ha line emission a 
5.2nm wide filter, centered at 657. Onm was used. For contin- 
uum subtraction we used a broad band red filter, centered at 
650nm. To cover the whole galaxy, a mosaic of 13 positions 
was obtained. The single maps were corrected for the image 
distortion of the focal reducer (with flux conservation) and then 
aligned using stellar DSS positions and merged to an Ha map 
of 55'x35' size. The flux calibration for the Ha map was done 
using the total line flux value given by Devereux et al. (1997) 
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Fig. 1. FIR maps of M33 (area 33'x33', North is up, East is left) at 60 (left top), 100 (right top), and 170yum (left bottom). The 
nOjjm map is restored to a spatial resolution of 100", while the resolutions at 60 and 100 ^/m are 60" and 80" respectively. For 
comparison, an R image of the galaxy is also shown. (The full resolution images of Fig. 1 can be obtained by anonymous ftp from 
ftp.mpia.de /pub/iso-arch). 



and using the Ha line fluxes of the prominent H n regions tab- 
ulated in Hippelein (1986), assuming a mean relative intensity 
of the [Nil] lines of 20% (McCall et al. 1985). 

3. Results 

In this section we will first present the general morphology of 
M 33, then we present the FIR photometry of discrete sources 
and their FIR spectral properties. 



3.1. Morphology 

Fig. 1 presents the maps obtained at the three FIR wavelengths, 
rebinned to a pixel size of 5". The spatial resolutions are 60" at 
60/mi and 80" at 100/im, as derived from the profile of IC133, 
the bright region in the NW edge. In the case of 170/mi the 
spatial resolution of ^120", as measured from the profile of 
NGC604 (IC133 is too faint at 170yum), was improved to 100" 
by deconvolving the original map with a PSF of 120". 
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Table 1. Comparison with IRAS results of Rice et al. (1990). 
<D are the aperture diameters used for the ISO photometry, Fgo 
stands for flux density F v at 60 //m in Jy, etc. 



Name 


O 






^60 


F ioo 






ISOPHOT 




IRAS 


total 




560 


1250 


420 


1260 


B651 


100" 


1.0 


1.9 


1.7 


4.5 


IC133 


140" 


6.3 


7.3 


11.8 


10.7 


IK60 


120" 


2.4 


4.6 


2.4 




NGC604 


260" 


29.5 


46.4 


33.6 


54.5 


IC131 


140" 


3.2 


6.4 


2.7 


3.2: 


NGC595 


140" 


13.0 


22.9 


12.8 




B220/1 


100" 


2.4 


5.4 


1.9 




B255/7 


100" 


1.7 


4.1 


1.6 




B248-51 


120" 


2.6 


5.5 


2.7 


8.3 



Fig. 2. Profiles through the center of M33 along position angle 
115° for the three FIR wavelengths. Here, the 60/im map is 
convolved with a Gaussian of 60" FWHM to achieve a similar 
spatial resolution as the 100/mi (80" FWHM), and the 170//m 
profile is taken from the Lucy-Richardson restored map. 



All three FIR maps clearly show the spiral arm structure of 
the galaxy and a large number of distinct sources. In addition 
to the spiral structure, an extended underlying component can 
be recognized, which becomes more significant with increas- 
ing wavelength. This is illustrated by the intensity profiles in 
Fig. 2, showing cuts through the three FIR maps at a position 
angle of 115°, along the minor axis. The structure of the pro- 
files is identical for the three filters, except that at the edge of 
the galaxy the 170/mi profile seems to decay slower than the 
others. 

The total flux density of M33 is 560Jy, 1250Jy, and 2200Jy 
at 60, 100, and 170/mi, respectively. 

3.2. FIR photometry of discrete sources 

About sixty sources were isolated from the FIR maps shown 
in Fig. 1. This number is large enough to allow statistical stud- 
ies for the FIR properties of extragalactic star-forming regions. 
Aperture photometry was done for these sources, with aperture 
diameters between 100" for small regions, up to 260" for the 
giant star-forming region NGC 604. Since in the crowded ar- 
eas the diffuse emission from M 33 cannot be unambigously 
determined, and since the introduction of a model for this com- 
ponent would be rather subjective and lead to arbitrary results, 
the photometry is performed without a subtraction of the dif- 
fuse emission. 

From IRAS maps of M33, Rice et al. (1990) localized 
a number of H n complexes and derived flux densities at 60 
and 100/im. Their 60/im flux densities (Table 1) are in reason- 
able agreement with ours, while their 100/im data are generally 
higher, probably due to the inferior spatial resolution of IRAS. 



3.3. Spectral energy distribution 

Fig. 3 depicts the SEDs as derived from the photometry for a 
number of selected areas. In the cases of the total galaxy and of 
NGC 604 and NGC 595 the IRAS (Rice et al. 1990) points are 
overplotted. The COBE/DIRBE (Odenwald et al. 1998) data 
overplotted in the SED for the total galaxy are somewhat lower. 
The explanation for this difference is that M33 is situated in a 
small hole in the cirrus foreground, and the off-galaxy back- 
ground measured by COBE/DIRBE with its large beam in the 
neighbourhood is probably too high, leading to lower flux den- 
sities for the galaxy. 

In all cases (Fig. 3) it is obvious from the distribution of the 
three FIR data points that the SEDs cannot be represented by 
a single blackbody, but need a combination of two blackbody 
emission components. We associate these two emissions with 
the two main morphological components measured directly in 
this paper, namely the localised (mainly warm) and the diffuse 
(mainly cold) component. Two blackbodies have four degrees 
of freedom and as there are only three data points, no unique 
combination can be derived. In the cases of M 33(total) and 
NGC 604 however, the IRAS data point at 25/«n puts an upper 
limit to the temperature of the warm blackbody component. 

At the long wavelength end Chini et al. (1986) found for 
Sb and Sc type galaxies observed by IRAS an average ratio 
Si200//m/S 60/jm ~ 0.034, with a scatter of a factor 2.5. Adopting 
the same ratio for M33 we insert a 1200 //m data point into the 
M33 SED. 

Attributing 30% of the signal at 25 /mi to hot dust, the fit 
with two modified blackbody functions (assuming an emissiv- 
ity of /3 oz AT 1 ) to the four data points between 25 and 170 /an 
yields T warm = 46.0 + 0.8 K, and T co , d = 16.7 + 0.5 K. As 
seen in the plot the cold blackbody curve comes very close to 
the 1200 /mi point estimated above. Adopting the temperature 
of the warm component to be rather constant over the galaxy, 
we keep T warm = 46 K fixed for fitting the cold component 
in the other areas. The resulting T co w varies only slightly, tak- 
ing values between 15.3 and 17.1 K, as demonstrated in the 
figure, with error bars of typically +0.5 K. This indicates that 
(coincidentially) the temperature of the diffuse cold dust com- 
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Fig. 3. Spectral energy distribution for total M33 (top left), for the central area (top center), for two prominent Hn regions 
(NGC604, NGC595), and for two areas in the interarm regime, W int. (offset 6' W, 5' S from the nucleus) and int. E (8' E 
from the nucleus). Dots are ISOPHOT data, circles stand for IRAS, asterisks for COBE/DIRBE data, and the plus with the large 
error bar at 1200/im is derived from flux density ratios as estimated by Chini et al., 1986). Two-component blackbody fits for 
temperatures of T warm = 46.0 K, T co u ~ 16.5 K, respectively, and with emissivity ft oc A~ 2 are overlaid as dashed and solid lines. 
For comparison, the dotted line shows a fit for/3 oc A 1 with T warm — 52.0 K, T co /j = 19 K. 



ponent is rather close to the temperature of the localised cold 
dust component. If one does the same for/? oc A -1 , one would 
get T warm = 52.5 ± 1.0 K, and T coM = 19.4 ± 0.7 K for M33 total , 
with a submm value a factor of 3 above the adopted one (dotted 
curve), while the variation of T coW was again rather small over 
the field. 

Therefore, we interpret the spatial variation of the flux den- 
sity ratio F„ arm /F co id as being due to a variation of their ampli- 
tudes, rather than being due to a change of their temperatures, 
which seem to remain almost constant over the galaxy. Since 
the warm and the cold component peak very close to 60 and 
170yum, respectively, the ISOPHOT 60 and 170yum flux densi- 
ties can be used as tracers of the warm and the cold dust. 

4. Discussion 

4.1. Spatial distribution of the warm dust 

In Fig. 4, the ratios F^q / Fno for all FIR knots are plotted versus 
the deprojected radial distance from the center of the galaxy. 
For the deprojection an inclination of 56° and a position an- 
gle of PA= 23° was assumed for the galaxy (Tully 1988). In 
addition a number of positions in the interarm regions (open 
circles) were selected. Firstly, it is obvious that the flux den- 
sity ratio for the interarm areas is lower than that of the bright 
regions at the same radius. Secondly, both the ratios in the in- 



terarm and bright regions decrease with increasing radial dis- 
tance, except for the bright H n regions NGC604 and 1013 3. In 
the case of the interarm regions, this decrease can be explained 
by a decrease in the 60/170 micron colour temperature of the 
diffuse emission of M 33. In the case of the bright star-forming 
regions, the decrease in the F 60 /F 110 ratio is attributed to an 
increasing relative contribution of the diffuse 170/im emission 
from the disk in the fainter, outer parts. Very bright regions, 
such as NGC604, but also those in the central area are almost 
unaffected by this diffuse disk component. A special case is 
IC133 due to its extraordinarily high Feo value. Schulman & 
Bregman (1995) found it associated with a bright X-ray source 
located in a hole in the Hi layer of the galaxy, indicating en- 
ergetic stellar winds and supernovae from massive stars. The 
average F(,q/F 170 ratio for the 20 bright star-forming regions 
is 0.32, with a dispersion of 0.05, and 0.165 for the interarm 
regions with a scatter of 0.03. 

The observed FIR brightness distribution is a superposi- 
tion of localized emission from star-forming regions and emis- 
sion from a diffuse component. The localized FIR emission 
is mainly powered by UV photons, whereas the diffuse FIR 
emission is powered by both UV and optical photons. Models 
of the UV-optical/FIR-submm SEDs which self consistently 
calculate a continuous distribution of dust temperatures based 
on radiative-transfer calculations (see Popescu & Tuffs 2002a) 
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predict that, for A < lOOyum, the FIR emission is dominated 
by the localized component. These models also predict that at 
longer wavelengths most emission originates from the diffuse 
disk component, and that both the localized and the diffuse dust 
emissions are predominantly powered by the UV light. 

We will try to separate the localized FIR emission (from 
star-forming regions) from the diffuse FIR emission compo- 
nent. While for central distances larger than r = 8' the interarm 
areas are wide enough to allow an estimate of the diffuse radi- 
ation from the disk, this is not possible in the inner area of the 
galaxy. 

Two simple approaches are performed. Firstly, assuming 
an exponential disk with a central flux density of 25 MJy/sr at 
60//m and a scale length of 6' (1.44kpc) for the diffuse emis- 
sion component would yield a residual map for the localized 
emission with a level in the inter-arm regions close to zero. 
The contribution of the diffuse disk to the total 60/im radiation 
would be 42%. This number agrees with the model for the Sb 
galaxy NGC 891 by Popescu et al. (2000) for which the con- 
tribution of the diffuse emission to the total is about 40% at 
60 //m. 

A second approach to obtain the contribution of the star- 
forming regions at 60 microns is to make a linear combination 
of the 60 and 170 /im maps. If we denote the diffuse emission 
component at 60 and 170 micron by and Ff JQ , respectively, 
and correspondingly the localized emission components by F l 6Q 
and F l UQ , then the arm and interarm colour ratios obtained pre- 
viously can be written as /' = F l 6f) IF l m and f d = F d Q /F d 70 , 
respectively. In this case the localized emission component at 
60 micron is given by the following equation, with F(,o and Fno 
in Jy: 



'60 



■fxF 



170 



60 



1- f/f 



(1) 



For f d = 0.165 and /' = 0.32, we find that 80% of the emis- 
sion at 60 micron comes from the localized component. So 
only 20% of the emission is diffuse at this wavelength. The 
20% value is lower than that obtained with the first method. 
However this is to be expected, since here we neglected the 
radial dependence of the f d and /' factors. The result in Fig. 

5 shows the star-forming regions well isolated against a zero- 
level background. 

4.2. Correlation of localized 60 /urn emission with Ha 
emission 

The location of the FIR knots is also traced in the distribution 
of H ii regions, best seen in the optical, in the light of the Ha 
line. In order to correlate the dust emission heated by the UV 
radiation from star-forming regions with the Ha emission from 
the associated H n regions, we make use of the map in Fig. 6. 
Before, the Ha map needs to be convolved with a Gaussian pro- 
file of 60" FWHM to achieve a similar spatial resolution. Fig. 

6 shows the central 33'x33' area of this map. Its morphology is 
almost indistinguishable from the map of localised warm dust 
emission at 60/j.m, F' 6Q (Fig. 5), except for the fact that the ratio 
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Fig. 4. The F^/ Fno ratio versus deprojected distance from the 
galaxy center. Large dots are for star-forming regions with high 
100//m flux density, small dots for faint ones. Some promi- 
nent regions are labeled by their names. Empty circles stand 
for inter-arm regions. 



V 




Fig. 5. Distribution of the localised warm dust component at 
60yum, F' 6Q , visualised by the scaled difference map 2(F&) - 
0.165 x Fno), with the factor 0.165 given by the average flux 
density ratio F(,o/Fno in the inter-arm regions, and the factor 2 
given by Eq. 1 . 



of FIR to H a emission appears to be higher in the central area 
than in the outskirts. 

The Fua/F 1 ^ ratio (Fig. 7, where F Ha is the flux in 
the line) for the FIR knots shows a clear systematic in- 
crease with increasing radial distance from the center (allow- 
ing for the [Nil] line contribution decreasing with distance, 
the slope would be even steeper). The low data point at r = 
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Fig. 6. Ha map of M33 convolved to a resolution of 60" 
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Fig. 7. FHa/F l m ratio for star-forming regions versus distance 
from the galaxy center radius. Symbol sizes indicate the 
brightness of the sources. (The ordinate is scaled in such a way 
that the mean ratio is 1) 



15 arcmin in Fig. 7 originates again from the star-forming re- 
gion IC133 with its unusually high 60yum flux density value. 
The least squares fit to the data points provides a slope of 
d(log(F Ha /^ ))/d(r) = 0.038/arcmin. 

A significant increase of the Fua/FpiR ratio towards the 
outer parts of the galaxy was also reported by Devereux et al. 
(1997). Their slope is however less steep, since it includes not 
only the star-forming regions but also locations inbetween the 
spiral arms, where the ratio is systematically lower. 



4.3. Extinction 

While the foreground extinction towards M 33 is negligible (A v 
= 0.09mag, de Vaucouleurs et al. 1991, or A v = 0.18mag, 
Burstein & Heiles 1982 and Schlegel et al. 1998), Israel & 
Kennicutt (1980) and Berkhuijsen (1983) derived from a com- 
parison of Ha data with radio flux densities at 21 and 6.2 cm 
extinction values of up to 2 mag. They also found a clear de- 
crease of the extinction for the Hn regions of M33 with in- 
creasing distance from the galaxy center of the form of Ay = 
2.3 - 0.07 r mag (with r in arcmin). 

From a pixelwise comparison of the Ha emission with a 
6 cm radio map Devereux et al. (1997) found no such trend of 
the internal extinction with radius but a scatter around a typi- 
cal value of the order of A v = 1. They therefore assigned the 
increase of the F Ha /F FIR ratio to a metallicity gradient in the 
sense, that hotter and more massive stars occur preferentially 
in the outer, metal-poorer regions of the galaxy. It seems that 
the radius dependence holds for the H n regions, but not for 
the diffuse emission, since, by means of multiwavelength ob- 
servations of Balmer and Paschen emission lines in several H n 
regions, Petersen & Gammelgaard (1997) again derived extinc- 
tion values decreasing with distance from the galaxy center in 
a similar way as found by Israel & Kennicutt (1980) and by 
Berkhuijsen (1983). 

When correcting the Ha fluxes according to the above Ay 
- r relation, the radial increase in Fiia/F 1 ^ mostly disappears, 
suggesting that this increase is mainly due to varying extinc- 
tion. This suggestion is supported by the relation between ex- 
tinction and cold dust radiation discussed in the next section. 

There might be other processes contributing to the trend 
shown in Fig. 7, such as a dust depletion at locations of low gas 
density, or/and effects due to metallicity gradients (Devereux et 
al. 1997). 



4.4. The cold dust 

Analogous to the isolation of the localized emission component 
at 60 micron (Sect. 4.1), we now construct a map of the diffuse 
cold dust component by suppressing from the 170/im map the 
contribution from the star-forming complexes: 



170 



^170 ~ f 60// 
l-fd/f 



(2) 



The scaled difference map is shown in Fig. 8. 

Since the F 60 /F ll0 ratio is higher than average in the bright 
knots (see Fig. 4), the difference map is overcorrected at these 
locations, thus yielding holes. Apart from these defects one can 
see that the spiral structure has almost disappeared. Instead, a 
disk like feature with an E-W extension of ~ 25' becomes ap- 
parent, revealing the distribution of the diffuse cold dust emis- 
sion. The 170jum-bright spiral arm segment in the South co- 
incides spatially with a region of very strong CO emission 
(Lequeux, priv.comm.), probably a region of forthcoming vi- 
olent star formation. 

The original 170//m map shows that large amounts of dust 
are located in or around the star-forming regions. This cold dust 
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Fig. 8. Distribution of the diffuse cold dust component at 
170 yum, Fj 70 , given by the scaled difference map 2(Fi7o - 
F 6 o/0.32). 



associated with the massive molecular cloud complexes is ex- 
pected to weaken the optical line emission of the embedded H n 
regions. A comparison between the Fua/F 1 ^ and the 170/im 
flux density shows (with the exception of NGC 604) that in- 
deed the brighter 170/luii sources exhibit lower F Ha /F6o ratios. 
Between the center of M33 and a deprojected radial distance of 
15' the F Ha /F l 6Q ratio increases by a factor of about 3.3, while 
the surface brightness S no (averaged over annuli) decreases by 
about 120MJysr _1 over the same distance. Assuming half of 
the cold dust to be located in front of the H n regions, and taking 
this dust to be responsible for the reddening and the variation 
of Fua/F l 6Q , this gradient would imply 

£ B -v ~ 0.01 x S no [MJy sr" 1 ] (3) 

Schlegel et al. (1998) derived from galactic cirrus extinc- 
tion a similar correlation between E(B-V) and the cold dust 
emission, £b-v = 0.018 x DT , where DT represents the 
IRAS 100/im flux density corrected to a reference temperature 
of 18.2K using the DIRBE/IRAS maps. 

The supposition that the intrinsic extinction is traced by the 
cold dust would give a straightforward explanation for the ap- 
parent contradiction between Devereux et al. (1997) who find 
no gradient in their pixel-wise comparison of F Ha and F^ cm 
(their figure 7) and the gradient for the H n regions presented 
by Israel & Kennicutt (1980), since the cold dust is concen- 
trated towards the star-forming regions as well. 

The cold dust constitutes a large amount of mass in galax- 
ies (e.g. Klaas et al. 2001, Popescu et al. 2002). Assuming an 
average grain size of 0. I/mi and a density of 3 g cirT 3 , the cold 
dust emission in M33 is equivalent to a mass of 6 x 1O 6 M . 
This is about 20% of what is found in M 3 1 (Haas et al. 1998), 
which is expected considering the higher mass of the latter. 



4.5. Star formation 

In the following paragraphs we will discuss how far the FIR 
and optical star formation calibrators lead to consistent results 
for the whole galaxy as well as for the individual star-forming 
regions. 

An integration of the two blackbody fit for M33 (see Fig. 
3) over the wavelength range 40 to 250//m, thus including the 
strong diffuse cold dust component, leads to 7.6 10~ n WirT 2 
= 1.5 10 9 L for M33. This is almost two times the value that 
one would derive from our 60 and 100//m points using Helou 
et al.'s (1988) relation F40-120 = 1-26 (2.58F 60 + F m ) 10~ 14 
[W m -2 ], with Feo and F100 in Jy, which do not include the cold 
dust component. 

Kennicutt (1998) derived from the star-burst synthesis 
models of Leitherer & Heckman (1995), for solar abundances, 
Salpeter IMF, and assuming that the dust re-radiates all the 
bolometric luminosity, a relation between star-forming rate and 
Lfir of 1.7 10~ 10 L F ir M yr~' L" 1 . Since the percentage of re- 
radiated light for normal galaxies is only around 30% (Popescu 
& Tuffs 2002b), the calibration factor used for starburst galax- 
ies will probably underestimate the SFR in M33. Hughes et al. 
(1998) derived for high redshifted HDF galaxies a calibration 
factor of 2.3 10~ 10 L FIR M yr~' L ', resulting in a star-forming 
rate of 0.35 M yr~' for M33, while the calibration of Rowan- 
Robinson et al. (1997) of 4.5 10 -10 L F1R Moyr 1 L" 1 yields a 
value of 0.69 M yr _1 , and a recent empirical relation between 
L FIR and star-forming rate (Bell 2003) gives 0.33 M yr _1 . 

Applying the Fn a -SFR calibration given by Kennicutt 
(1998) to the total Ha line luminosity of M 33, and correcting 
it for a mean A v = 1.2, (Israel & Kennicutt 1980, Devereux et 
al. 1997) and for a mean [N 11] contribution (averaged over the 
galaxy) of 20% (McCall et al. 1985, Vilchez et al. 1988), the 
star-forming rate results to S FR(Ha) = 0.47 M yr _1 , which is 
in satisfactory agreement with the FIR values. 

When we, however, pick a substructure of the galaxy - for 
instance NGC 604, the brightest Hn region in M 33 - no con- 
sistency can be achieved between the SFRs derived from the 
optical and from the FIR. For this example, the Ha luminos- 
ity corrected for an intrinsic extinction of Ay = 0.8 (the mean 
of Israel & Kennicutt 1980, and of Peterson & Gammelgaard 
1997), and for 17% [N 11] contribution for NGC 604 (Vilchez et 
al. 1988), yields a star-forming rate of 0.046 M yr~'. By com- 
parison, if all the UV light were locally absorbed, its observed 
FIR luminosity of 5.8 1O 7 L would yield SFR values of only 
between 0.010 (Kennicutt 1998) and 0.026 M yr" 1 (Rowan- 
Robinson et al. 1997), just 40% of the value derived from the 
optical. Thus, in this case, some 60% of the stellar UV light 
escapes the dust sphere around the H 11 region and contributes 
to the heating of the diffuse dust in the galaxy disk! 

4.6. A local radio-FIR correlation 

One of the tightest correlations in extragalactic astronomy is 
that between the FIR luminosity, as calculated from the IRAS 
flux densities at 60 and 100 /im (for example combined in the 
above-mentioned relation F40-120 by Helou et al. 1988), and 
the differential radio continuum luminosity in the GHz range 
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for late-type galaxies. It was originally found by de Jong et al. 
(1985), Helou et al. (1985) and Wunderlich et al. (1987), and 
holds for the spatially integrated emissions of objects whose 
luminosity is dominated by star formation and not by an active 
nucleus. 

This radio-FIR correlation can be theoretically explained, 
since the Interstellar Medium in such galaxies acts as a global 
calorimeter for the relativistic electrons and, in large part, for 
the U V photons, both of which are generated by the stars ( Volk 
1989). Usually the correlation is found to be somewhat steeper 
than unity, the radio luminosity increasing faster than linearly 
with the FIR luminosity. Restricting the FIR to the IRAS 60 fim 
luminosity, Yun et al. (2001) have recently obtained a nearly 
linear correlation, based on a very large sample of 1809 galax- 
ies from the local Universe, which Bell (2003) considers as a 
matter of conspiracy of different conflicting effects. 

Our present work allows us to investigate whether such a 
correlation also exists for the typically sub-kpc size galactic 
substructures such as the star-forming regions in M33. For this 
purpose we made use of the 6.3 cm radio map of Buczilowski 
& Beck (1987), which has a spatial resolution of 2.4'. Since 
this resolution is considerably lower than the resolution of our 
60 and lOOyum maps we smoothed the FIR maps to a compa- 
rable spatial resolution, before rederiving the photometry for 
the knots in M33. For the estimate of the FIR fluxes we used 
the relation given by Helou et al. (1988) to allow a comparison 
with IRAS data. Fig. 9 shows the radio flux plotted versus the 
FIR flux. There appears to be indeed a correlation. Even taking 
into account the difficulties inherent in the data reduction, the 
slope of this correlation is smaller than unity, of the order of 
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Fig. 9. Plot of the monochromatic radio luminosity versus the 
FIR luminosities for M33 and its star-forming regions (filled 
circles) together with the data for the Effelsberg 100-m galaxy 
sample (Wunderlich et al. 1987, open circles). The dotted line 
has a slope of 1.10 (Wunderlich & Klein 1988). 



0.9. We can ask ourselves why we should expect such a "local" 
correlation. 

A priori this is not obvious. The absorption of the non- 
ionizing stellar UV photons depends strongly on the dust col- 
umn density through the star-forming regions. And even for 
NGC 604 about 60 percent of the UV light escapes the dust 
sphere of the region, as mentioned above. For smaller regions 
this fraction is expected to be higher. The rest goes into the 
diffuse dust in the galactic disk. This is the dominant FIR part 
in the standard global radio-FIR correlation, but it can be ne- 
glected in the "local" correlation, unless the surface filling fac- 
tor of the ensemble of regions becomes a significant fraction, 
say more than 10 percent of the disk area. Nevertheless, the 
fact that we did not include the 170/luii in the FIR luminosity 
of Fig. 9 emphasizes the shorter wavelengths, and we noted 
earlier that a large fraction (about 80%) of the emission at 
60 ^m comes from the localized component. Thus, some frac- 
tion, say 30%, of the spatially integrated IRAS FIR luminos- 
ity may come from individual star-forming regions and will be 
proportional to their respective star formation rates. In fact, a 
higher star formation rate will probably imply a higher dust 
column density and thus a stronger than linear dependence of 
the FIR luminosity on the local star formation rate. 

Let us next consider the radio emission from individual 
regions. If we assume that the sources of the radio emitting 
relativistic electrons are mainly the remnants of core collapse 
supernovae, then their progenitor stars will have remained in 
the star-forming region during their evolution. In addition to 
this non-thermal radiation, the progenitor stars will give rise 
to local thermal (free-free) radio emission. Therefore the nec- 
essary condition for a "local" correlation is satisfied, because 
the sources for both emissions are localized in the star-forming 
region. 

Also, for a face-on galaxy, the star-forming regions are 
detected along the direction perpendicular to the disk. Along 
that direction we see the integrated emission, and the elec- 
trons confined to this column will have lost their energy by 
Inverse Compton and synchrotron emission, and thus their lu- 
minosity can be proportional to the local star formation rate. 
Much more important is the following effect: the supernova 
remnants which accelerate the electrons, confine them during 
the active remnant life time of ~ 10 5 yr, during which time 
they radiate with a very hard energy spectrum. Acceleration 
automatically implies confinement, because in the shock frame 
the particles are convected in the downstream direction. During 
their confinement the electrons radiate in the interior effective 
magnetic field which can be expected to be several times, say 

5 times higher than in the average interstellar medium (Bell 

6 Lucek 2001; Berezhko et al. 2002). In agreement with ob- 
servations (Biermann 1976, Ulvestad 1982) this implies a ra- 
tio of ~ (25 x 10 5 yr)/(3 x 10 7 yr) ~ 10" 1 between the syn- 
chrotron luminosity from the ensemble of star-forming regions 
to the diffuse luminosity from the overall galaxy where the par- 
ticles have a life time of ~ 3 x 10 7 yrs. Nevertheless, most 
of the relativistic electrons will leave the region after their 
release from the sources before having lost their energy ra- 
diatively. This is easily seen, since their radiative loss length 
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(k(E) x ti oss (E)) 1/2 > few kpc 1 exceeds the typical sizes of 
star-forming regions, and this is even a lower limit (Ptuskin et 
al. 1997). Therefore, like the cold dust emission from Fig. 8, 
the diffuse synchrotron emission will constitute a featureless 
disk-like structure and can be neglected in the local correla- 
tion. In fact, the escape from a star-forming region will be the 
more effective the higher the cold gas mass is there, and should 
therefore decrease with increasing FIR luminosity, if anything. 

Therefore, we must expect some local radio-FIR correla- 
tion in star-forming regions, and even one which has a slope 
< 1. This explains the observational results shown in Fig. 9. 
However, this correlation is primarily due to the confinement 
of synchrotron electrons in their sources (supplemented by lo- 
cal free-free emission), and due to local dust absorption. As 
such it will primarily be seen through the warm dust emission 
observed by IRAS. For the total FIR as measured by ISOPHOT 
we expect a much weaker correlation, if at all. Quantitatively, 
the correlation depends very much on the local transport con- 
ditions for photons and relativistic electrons and has little con- 
nection with the global radio-FIR correlation because the in- 
gredients of the latter correlation can be neglected in the lo- 
cal correlation. In this picture one would predict that the radio 
emission responsible for the local correlation will have a hard 
spectrum with a spectral index flatter than about 0.5. 

5. Conclusions. 

The FIR radiation from M 33 is composed of a warm (~ 45K) 
and a cold (~ 17K) dust emission component, which were sep- 
arated by fitting the SEDs of the entire galaxy and of individ- 
ual regions with two modified Planck functions of emissivity 
P cc A' 2 . We spatially isolated the two dust emission compo- 
nents and proved that they can be identified with separate mor- 
phological components. The warm dust is preferentially asso- 
ciated with the star-forming regions. The cold dust arrises in 
part from a diffuse disk of about 25' diameter, but also has a 
localized component associated with the molecular clouds in 
star-forming regions. 

Estimates of the star formation rate derived from the FIR 
for the whole galaxy are broadly consistent with that derived 
from the Ha emission. In the case of individual star-forming 
regions, a consistency between rates derived from the optical 
and from the FIR allows only a fraction of the UV radiation to 
be absorbed locally. 

The correlation between the F Ha /F^o ratio and the 170 //m 
surface brightness is in agreement with the extinction relation 
established by IRAS/DIRBE observations for the cirrus in the 
Milky Way. 

The ISOPHOT maps were also used to derive a local radio- 
FIR correlation for star-forming regions in M 33. It is argued 
that this local correlation is due to quite different emission com- 
ponents than to those that lead to the well-known global radio- 
FIR correlation. 



1 Here k{E) is the spatial diffusion coefficient of electrons of kinetic 
energy E, and T/ .„(£) is the electron energy loss time. The product 
k(E) x t /o ,„(£) oc E 0A is almost energy-independent. 
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